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Introduction 


The electrical activity of muscle has long been studied by 
recording from the surface of a muscle and by recording the 
activity of individual motor units with concentric needle elec- 
trodes. Both techniques are simple and reliable enough to be 
used routinely In the kinesiology and in the diagnosis of many 
diseases tf the muscles and their motor nerves [l]. 

The Interference pattern electromyogram (EMG) which is ob- 
served when recording from the muacle*a surface is a summation of 
action potentials of a number of motor units. The information 
about these motor units latent in it Is not readily apparent 
because of the very complex pattern. Consequently, the surface EMG 
has not yet seen widespread clinical application. The purpose of 
this study is to investigate the quantitative regularities of 
Interference pattern formation by motor unit action potentials. 

To this end, we have considered the parameters of a single motor 
unit and how they relate to the Fourier transform analysis of an 
EMG. The Fourier transform of the simulated electromyogram is 
compared with the Fourier transform of the actual EMG recorded 
from various human muscles using surface electrodes. 

The Parameters of the Motor Unit 

The steady state discharge of single neurons is often charac- 
terized by a probability density function for the intervals 
between successive nerve impulses called the Inter Spike Intervals 
(ISI). 

Buchthal and his associates [2,3]* Clamann [4], Leifer [5] 
and several other investigators have shown that the ISI of suc- 
cessive motor unit spikes are normally distributed. The rates 
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of steady state discharge in the cats are usaully between 5 to 30 
pulses per second (pps) [6], 

Clamann [4] has reported the firing rates between 7 to 25 pps 
in the human brachial biceps. The transient discharge rates may 
be higher, pert icularly in smaller muscles. A peak frequency of 
150 pps has been reported in the adductor pollicis in the human 
hand [7]. In unfatigued human muscles, the mean and the standard 
deviation of the probability density function characterizing the 
131 are functionally related [4], making it possible to completely 
specify the stochastic process by determining its mean. For the 

human brachial biceps muscle he arrives ar the following functional 

2 

relation between the mean (u) and the variance ( a ) of the ISI: 
o » 9*1 x 10 ^ + 4.0 msec (1) 

When successive intervals are uncorrelated, the probability 
density function and its mean and variance completely specify the 
sta 1st leal process generating the intervals. In the normal, 
unfatigued muscle, under conditions of low or moderate activation, 
zhe steady discharges of motoneurons generally do not show any 
dependence between successive intervals [4,5*8 >9] • However, at 
high levels of muscular activity and particularly during fatigue, 
the discharges tend to be grouped so that successive intervals 
will nc longer be uncorrelated. 

The EMG of a single motor unit can be looked at as a time 

function defined uy a convolution integral: 
t 

e(t) » J h(t-x) p(t )dt (2) 

where p(t) is a point process (i.e., a series of unit impulses 
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or Dirac delta functions) which passes through a filter whose 
Impulse response ic h(t)&s in Fig. 1. Tb3 time function h(t) des- 
cribes the shape of a single motor unit action potential. The 
width of the motor unit action potential Is usually between 3 to 
20 msec. [2]. The conduction process of the action potential 
along the muscle fiber is such that successive muscle action 
potentials recorded from the same motor unit cannot overlap. 

They remain distinct events. A typical EMG of a synthetic motor 
unit is shown in Figure 2. The time origin is arbitrarily fixed 
so that the first motor unit starts at t « i msec. The Fourier 
transform (FT) of equation (2) is 

E(Jw) - H(Jw) P(Jw) (3) 

where H(Jw) is the FT of the impulse response and P(jw) is the FT 
of the random pulse train. 

Let • x l + t 2 + and - 0. Then the ith pulse 

occurs at time {i + t^). For a pulse train of N M pulses, 
p(t) may be written as 


N 

p(t) - z 6(t-e~t.) 
i-0 1 


(4) 


and its FT is 


P(Jw) 


n ) 
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N 
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(5) 


Without any loss of generality, we will from now on consider 
i - 0, i.e., the time origin Is defined at the occurrence of the 
first pulse. Tberofore, 
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P(Jw) 


N -jwt, 
L e 
i«0 
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where tj Is a sum of i normally distributed independent ISI. 
Ensemble Average °l m »> 

The (i+l)th pulse occurs at time which is a sum of i 
Independent random variables 

t A = + t 2 + . . . + -r^ 

p 

The mean (u) and variance (a ) of each ISI is assumed bo be 
the same. Thus the mean and variance of t^ are given by: 

E ( ) = E (t 1 + Tp +...+ T^) 

i 


- C E(t . ) = iu 
J=1 J 

V(t.) = V(t 7 + T. +...+ T 


i' 


i 

* z 
J-l 


V(Tj) 


lo 


rhe probability density function of t^ is given by 
1 


f(t x ) 


✓2*1 


exp 


-(t^lu) 2 

2io 2 
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rhe Fourier transform P(Jw) is a random function of i » 1,2, 
. . . , N. The ensemble average is given by 


FTjwT * E [P(Jw, t^, tp,*.*, t^)] 

N -Jwt. 

* E [ I e *] 

1 = 0 

N -Jwt. 

« l E(e 1 ) (10) 

1=0 
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The expected value of e is given by 

-Jwt* J* 


-Jwt, j* -Jwt, 

E(e *) = f e 1 


f(t i ) dt x 


00 

— — /■ 
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-Jwt x -Oj-Im)' 
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Let -ip • x. Then 


-Jwt, 5 -Jwip 
E[e 1 ] = 

/27T 


-jwin 


e =-5 — dt A 

2io z 1 
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cos wx e 2io dx 
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( 11 ) 


■ e *e 

The ensemble average P(Jw) from equations (10) and (11) is given 
by 

PTjwT -re (12) 

1-0 £ 


From equation (2) and (12), the ensemble average of the Fourier 
Transform of a single motor unit is given by 

eTJmT - FTJw) H(Jw) '.13) 

and 

|e7TwT I = |P(Jw) I | H( j w) | (14) 

In Figure 3> the FT of the EMG given by equation 14 is plotted for 
u = 20 msec, a ® 4.36 msec and N * 50 (which represents about 50 
pulses on the average in one second interval*. H(jw) is taken 
for a diphasic pulse with peak to peak height 2b and width 2c as 
shown in Figure 2. The FT for this pulse 1 b 

H(Jw) - ^5 e‘ Jwc (2s In (*§) - sin (wc)) (15) 

cw 
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EMG Simulation 

The EMG simulation was carrieu out on an IBM 1800 process 
computer. Our simulation is very similar to the work done by 
Person and Libkind CIO, 111]. Their studies were, however, limited to 
the symmetrical biphasic motor unit potentials. 

In our model, M is the number of active motor units. " * 

irterspike intervals (ISI) are independent random values having 
& normal distribution with mean u and standard derivation a. 

Various shapes of motor unit potentials were used in this study 
with widths from 4 to 16 msecs. The pulse shapes were choser. to 
have a zero mean value. The ISI was always taken to be greater 
than the width of the pulse. 

The first pulse for the ith motor unit was taken at a random 
interval uniformly distributed between zero and u msec. All other 
I3 T were normally distributed. The interference pattern was 
arrived at by surarairg up the action potentials of active motor 
unit v' in the time interval from 0 to 1C 50 msec. The samples were 
1 msec apart. The Fourier transform was obtained using 1024 samples 
and employing a fast Pourier transform algorithm. 

Simulation Results 

With a suall number of motor units active, the pulses were 
frequently grouoed with a "pile up effect". With an increase In 
the number of active motor units the time Intervals free of 
impulses gradually diminished. A "saturation" level is reached 
with about 16 active motor units at mean frequency of 20 pps. 

Figure 4A shows that synthetic EMGs, especially the "saturated" 
ones, bears a close resemblance to the actual surface EMG, examples 
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of which appear in Figure 4B. 

The theoretical analysis indicates that the average FT is 
the product of the FT of the motor units pulse shape and the 
ensemble average of the FT of a randomly occurring impulse train. 
Figure 5 shows the | FT | for one motor unit with a symmetric 
biphasic pulse of 8 msec width as shown in Figure 2 Cp • 20 vec, 
a - 4.36 msec, M * 13. Figure 6 shows the average | PT | fcr 15 
such independent motor units. The smooth curve in Figure 6 is 
obtained from equation 14 for 30 spike sequence. The important 
point in this analysis is that average |FT| of several Identical 
motor units approaches the FT of a single motor unit's pulse 
shape. Figure 7 shows some samples of the average | FT | of the 
surface EMG data recorded from the gastrocnemius-soleus muscle. 
Discussion 

In this analysis it has been assumed that mean and variance 
of the interspike interval remains constant and also that the 
shape of the various motor unit potentials in a muscle are the 
same. The theoretical analysis shows that the avorage |FT| of 
the synthetic EMG approaches the | FT | of the motor unit potential. 
This has been confirmed by simulation studies except at the very 
low end of the spectrum. 

If one assumes that the mean and the variance of the 1ST for 
the real EMG remain nearly stable over a short interval during 
constant force and also that the shape of the various motor unit 
potentials within the recording field of the electrodes are the 
same, then It is theoretically feasible to calculate the shape of 
the motor u.*it potential from the average FT. The practicality 
of this procedure is now being investigated. 
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SYNTHETIC EMO 



FIGURE 1: EMO Model as low pass filter 
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FIGURE 2: Synthetic EMG for a single motor unit 

with the parameters u » 50 msec and a • 6.27 msec 
[Time « *00 msecl 


578 - 





o 250 500 


FREQUENCY (H x ) 


FIGURE 3: Theoretical absolute Fourier transform 

of the EMG for a symmetric blphasic pulse of 8 
m3ec width ,u * 20 msec, c « 4.36 msec and the num- 
ber cf pulses is 50. 



FIGURE 4A: Synthetic liMG generated for a symme- 

tric triphasic pulse. The four traces are for 
1, 5, 10 and 20 active motor units. 
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FIGURE 5: Absolute Fourier transform for a single 

active motor unit with biphat ic pulse shape (u * 

20 msec and a • 4.36 msec). 
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FIGURE 6: Average absolute Fourier transform for 

fifteen independent active motor units with bi- 
phaslc pulse shape ( v * 20 msec and o • 4,36 nsec). 
The ^mocth curve is the theoretical absolute 
Fourier transform. 
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FIGURE 7; Average absolute Fourier transform 
of surface EMC data from gastrocnemius-soleus 
muscle. The three parts are average of 1. 5 
and 15 samples. ' 


